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Abstract 

In this paper we study the class of differential operators T = X^=i QjD-' 
with polynomial coefficients Qj in one complex variable satisfying the 
condition deg Qj < j with equality for at least one j. We show that if 
degQfe < k then the root with the largest modulus of the nth degree 
eigenpolynomial p n of T tends to infinity when n — > oo, as opposed to the 
case when degQ^ = k, which we have treated previously in [2], Moreover 
we present an explicit conjecture and partial results on the growth of the 
largest modulus of the roots of p n . Based on this conjecture we deduce 
the algebraic equation satisfied by the Cauchy transform of the asymp- 
totic root measure of the appropriately scaled eigenpolynomials, for which 
the union of all roots is conjecturally contained in a compact set. 

1 Introduction 

In this paper we study asymptotic properties of roots in certain families of 
eigenpolynomials. Namely, consider a linear differential operator 

k 
3 = 1 

where D = d/dz and the Qj are complex polynomials in a single variable z 
satisfying the condition deg Qj < j for all j, and degQ^ < k for the leading 
term. Such operators will be referred to as degenerate exactly- solvable operators, 
see Definition 1 below. In this paper we study polynomial eigenfunctions of such 
operators, that is polynomials satisfying 

T(p n ) = X n p n (1) 
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for some value of the spectral parameter A„, where n is a positive integer and 
degp„ = n. 

The basic motivation for this study comes from two sources: 1) a classical 
question going back to S. Bochner, and 2) the generalized Bochner problem, 
which we describe below. 

1) In 1929 Bochner asked about the classification of differential equations JT]) 
having an infinite sequence of orthogonal polynomial solutions, see [TSj- Such 
a system of polynomials {pn}^ = Q which are both eigenpolynomials of some fi- 
nite order differential operator and orthogonal with respect to some suitable 
inner product, are referred to as Bochner- Krall orthogonal polynomial systems 
(BKS), and the corresponding operators are called Bochner-Krall operators. It 
is an open problem to classify all BKS - a complete classification is only known 
for Bochner-Krall operators of order k < 4, and the corresponding BKS are var- 
ious classical systems such as the Jacobi type, the Laguerre type, the Legendre 
type and the Bessel and Hermite polynomials, see [7]. 

2) The problem of a general classisfication of linear differential operators for 
which the eigenvalue problem (TT]) has a certain number of eigenfunctions in the 
form of a finite-order polynomial in some variables, is referred to as the gener- 
alized Bochner problem, see |19] and [20j . In the former paper a classification of 
operators possessing infinitely many finite-dimensional subspaces with a basis 
in polynomials is presented, and in the latter paper a general method has been 
formulated for generating eigenvalue problems for linear differential operators 
in one and several variables possessing polynomial solutions. 

Notice that for the operators considered here the sequence of eigenpolyno- 
mials is in general not an orthogonal system and it can therefore not be studied 
by means of the extensive theory known for such systems. 

Definition 1. We call a linear differential operator T of the fcth order exactly- 
solvable if it preserves the infinite flag Vo C V\ C Vi C • • • C V n C • • • , where 
V n is the linear space of all polynomials of degree less than or equal to rtQ Or, 
equivalently, the problem {T]) has an infinite sequence of polynomial eigenfunc- 
tions if and only if the operator T is exactly- solvable, see [21] . 

Notice that any exactly-solvable operator is of the form T = Y^j=i QjD 3 ■ 
They split into two major classes: non- degenerate and degenerate, where in the 
former case degQ^ = k, and in the latter case degQfc < k for the leading term. 
The major difference between these two classes is that in the non-degenerate 
case the union of all roots of all eigenpolynomials of T is contained in a compact 
set (see [5]), contrary to the degenerate case, which we prove in this paper. 

1 Correspondingly, a linear differential operator of the fcth order is called quasi- exactly- 
solvable if it preserves the space P n for some fixed n. 
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The importance of studying eigenpolynomials of exactly-solvable operators is 
among other things motivated by numerous examples coming from classical or- 
thogonal polynomials. Our study can be considered as a natural generalization 
of the behaviour of the maximal root for classical orthogonal polynomial families 
such as the Laguerre and Hermite polynomials, which appear as solutions to the 
eigenvalueproblem {T]) for certain choices on the polynomial coefficients Qj for 
a second-order degenerate exactly-solvable operator; the Laguerre polynomials 
appear as solutions to the differential equation zy"(z) + (1 — z)y'(z) +ny(z) = 0, 
and the Hermite polynomials are solutions to the differential equation y"{z) — 
2zy'(z) + 2ny(z) = where n is a nonnegative integer. Recent studies and 
interesting results on the asymptotic zero behaviour for these polynomials and 
the corresponding generalized polynomials can be found in e.g [6], [9].[l3]. [17], 
[18] , [11] , pj] and references therein. In [11] one can find bounds on the spacing 
of zeros of certain functions belonging to the Laguerre-Polya class satisfying a 
second order differential equation, and as a corollary new sharp inequalities on 
the extreme zeros of the Hermite, Laguerre and Jacobi polynomials are estab- 
lished. 

Let us briefly recall our previous results. In 2] we treated the asymp- 
totic zero distribution for polynomial families appearing as solutions to (fT]) 
where T is an arbitrary non-degenerate exactly-solvable operator. This seems 
to be a natural generalization to higher orders of the Gauss hypergeometric 
equation. As a special case, the Jacobi polynomials appear as solutions to 
(z 2 — l)y"(z) + (az + b)y'(z) + cy(z) = 0, where a, b and c are constants satisfy- 
ing a > b, a + b > and c = n(l — a — n) for some nonnegative integer n. It is a 
classical fact that the zeros of the Jacobi polynomials lie in the interval [—1,1] 
and that their density in this interval is proportional to 1/yT — \z\ 2 when the 
degree n tends to infinity, which follows from the general theory of orthogonal 
polynomial systems. However, for higher-order operators of this kind, the se- 
quence of eigenpolynomials is in general not an orthogonal system. In [2] we 
proved that when n — > oo, the roots of the nth degree eigenpolynomial p n for a 
non-degenerate exactly-solvable operator are distributed according to a certain 
probability measure which has compact support and which depends only on the 
leading polynomial Qk- Namely, 

Theorem A. Let Qk be a monic polynomial of degree k. Then there exists 
a unique probability measure fiQ k with compact support whose Cauchy trans- 
form C(z) = / dM z °_ fc c (C) satisfies C(z) k = 1/Q k (z) for almost all z E C. 

Theorem B. Let Qk and fiQ k be as in Theorem A. Then supp /iQ fc is the union 
of finitely many smooth curve segments, and each of these curves is mapped to 
a straight line by the locally defined mapping &(z) = f Qk(z)~ 1 / k dz. Moreover, 
supp /iQ fc contains all the zeros of Qk, is contained in the convex hull of the 
zeros of Qk and is connected and has connected complement. 
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If p n is a polynomial of degree n we construct the probability measure /i n 
by placing the point mass of size i at each zero of p n , and we call fi n the root 
measure of p n . The following is our main result from [5] : 

Theorem C. Let p n be the monic degree n eigenpolynomial of a non-degenerate 
exactly-solvable operator, and let fi n be the root measure of p n . Then fi n con- 
verges weakly to /iQ k when n — > oo. 

To illustrate, we show the zeros of some polynomial eigenfunctions for the 
non-degenerate exactly-solvable operator T = Q5D 5 , where Q5 = (z — 2+2i)(z+ 
l — 2i)(z + 3 + i)(z + 2i)(z — 2i — 2). In the pictures below large dots represent the 
zeros of Q§ and small dots represent the zeros of the eigenpolynomials pso,P75 
and piqq respectively. 



n=50 



n=75 



n=100. 



As a consequence of the above results we were able to prove a special case of 
a general conjecture describing the leading terms of all Bochner-Krall operators, 
see [5]. 

In the present paper we partially extend the above results to the case of de- 
generate exactly-solvable operators. Numerical evidence shows that the roots of 
the nth degree eigenpolynomial are distributed along a tree in this case too, but 
that the limiting root measure is compactly supported only after an appropraite 
scaling of the roots. We will assume wlog that p n is monic. 

We start with the following preliminary result: 

Lemma 1. Let T = ^j=i QjD-* be a degenerate exactly-solvable operator of 
order k. Then, for all sufficiently large integers n, there exists a unique constant 
X n and a unique monic polynomial p n of degree n which satisfy T(p n ) = \ n p n . 
If deg Qj = j for precisely one value j < k, then there exists a unique constant 
X n and a unique monic polynomial p n of degree n which satisfy T(p n ) = \ n p n 
for every integer n — 1,2,.... 

In what follows we denote by r n the largest modulus of all roots of the unique 
and monic nth degree eigenpolynomial p n of T, i.e. r n = max{|z| : p n {z) = 0}. 
These are our main results: 
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Theorem Let T be a degenerate exactly-solvable operator of order k. Then 
r n — > oo when n — > oo. 

Next we establish a lower bound for r n when n — > oo. 

Theorem 2. Let T = V*' , Q.IP = E*=i (EiS* ' a o,i zi ) Dj be a de 9 en ~ 
erate exactly-solvable operator of order k. Then for any 7 < b we have 

r„ 

lim — = 00, 

where 

+ 
min 



je[i,fc-i] \k - j + degQj ■- degQ fc / ' 

anrf where the notation min + means that the minimum is taken only over posi- 
tive terms (k — j + deg Qj — deg Qk)- 

Corollary 1. Let T = Ej=i QjD-* be a degenerate exactly-solvable operator 
of order k such that deg Qj < jo for all j > jo, and in particular degQfe = jo, 
where jo is the largest j such that degQj = j. Then linin^oo ^ = 00 for any 
7 < 1. 

Corollary 2. Let T = Ej=i QjD 3 be a degenerate exactly- solvable operator 
of order k such that deg Qj = for all j > jo, where jo is the largest j for which 
degQj = j. Then Hindoo ^ = 00 for any 7 < ^j^-. 

In fact our extensive numerical experiments and natural heuristic arguments 
(see Section 3) support the following 

Main Conjecture. Let T = Ej— 1 Qj^ be a degenerate exactly-solvable op- 
erator of order k and denote by jo the largest j for which deg Qj = j . Then 

,. r n 

lim — = ct, 

n — >oc Tl 

where ct > is a positive constant and 

j - jo 



d:= 



max 



jebo+i,fe] \j - degQ-j 
Remark. Note that Main Conjecture implies Theorem 2 since b < d. 
The next two theorems support the above conjecture: 



2 This theorem is joint work with H. Rullgard. 
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Theorem 3. Let T be a degenerate exactly-solvable operator of order k consist- 
ing of precisely two terms: T = Qj D^° + QkD k . Then there exists a positive 
constant c such that 

lim inf — ^ > c 

where d := max ie[io+a] ( j^gfej ) = fc-degQ, • 

This result can be generalized to operators consisting of any number of terms, 
but with certain conditions on the degree of the polynomial coefficients Qj for 
which j > jo 7 where jo is the largest j for which deg Qj = j. Namely, 

Theorem 4. LetT be a degenerate exactly- solvable operator of order k. Denote 
by jo the largest j such that deg Qj — j and let (j — deg Qj) > (k — deg Qk) for 
every j > jo ■ Then there exists a positive constant c > such that 

lim inf — ^ > c 

n^oo n 

where d := max j6[joW , t] ( , J ^ Q ) = k ^ Qk ■ 



Numerical experiments show that roots of eigenpolynomials scaled according 
to the Main Conjecture fill certain interesting curves in C. To illustrate this phe- 
nomenon let us present three typical pictures. Below p n denotes the nth degree 
unique and monic eigenpolynomial of the given operator, and q n (z) = p n (n d z) 
denotes the corresponding appropriately scaled polynomial, where d is as in 
Main Conjecture, and for which the union of all roots is (conjecturally) con- 
tained in a compact set. 

Fig.l: Fig. 2: Fig. 3: 



roots of roots of roots of 

9100(2) = _pioo(100z) qioaiz) = pioo(lOOz) £100(2) = Pioo(lOOz) 

Fig.l: 7\ = zD + zD 2 + zD 3 + zD 4 + zD 5 . 

Fig. 2: T 2 = z 2 D 2 + D 7 . 

Fig. 3: T 3 = z 3 D 3 + z 2 D A + zD 5 . 

In Section 3 we will derive the (conjectural) algebraic equation satisfied by 
the Cauchy transform C(z) of the asymptotic root measure of the scaled eigen- 
polynomial q n for an arbitrary degenerate exactly-solvable operator. From this 



G 



equation one can obtain detailed information about the above curves and also 
conclude which terms of the operator that are relevant for the asymptotic zero 
distribution of its eigenpolynomials@ Numerical evidence clearly illustrates that 
distinct operators whose scaled eigenpolynomials satisfy the same Cauchy trans- 
form equation when n — > oo, will yield identical asymptotic zero distributions. 
Below we show one such example. For further details see Section 4.3. 



T 4 = z 3 D 3 + z 2 D 5 , T 4 = z 2 D 2 + z 3 D 3 + zD A + z 2 D 5 + D 6 , 

roots of 3100(2) = pioo(100 2/3 z) roots of q W o{z) = pioo(100 2/3 z) 

Let us finally mention some possible applications of our results and direc- 
tions for further reasearch. As was mentioned earlier, operators of the type we 
consider occur in the theory of Bochner-Krall orthogonal systems. A great deal 
is known about the asymptotic zero distribution of orthogonal polynomials. By 
comparing such known results with results on the asymptotic zero distribution 
of eigenpolynomials of degenerate exactly-solvable operators we believe it will 
be possible to gain new insight into the nature of BKS. 

This paper is structured as follows. In Section 2 we give the proofs of the 
lemma, the theorems and the corollaries stated in this section. In Section 3 we 
explain how we arrived at Main Conjecture and how we obtained as its corollary 
the algebraic Cauchy transform equation. In Section 4 we give numerical evi- 
dence supporting Main Conjecture and its corollary. In Section 5 (Appendix) 
we give the detailed calculations led to the corollary of Main Conjecture, and we 
also prove that for a class of operators of the type we consider, the conjectured 
upper bound for r„ implies the conjectured lower bound. Finally in Section 6 
we discuss some open problems and directions for further research. 

Acknowledgements. I am sincerely greatful to my PhD advisor Professor 
B. Shapiro for introducing me to this very fascinating problem and for his con- 
stant support during my work. I would also like to thank Professor J-E. Bjork 
for stimulating discussions on the topic. I am also greatly obliged to anonymous 
referees for their interesting suggestions which contributed to the improvement 
of this paper. My research was supported by Stockholm University. 

3 As was mentioned earlier, in the non-degenerate case which we have treated previously, the 
asymptotic zero distribution of the eigenpolynomials depends only on the leading coefficient 
Qj,. For the operators considered here however, the situation is more complicated. 
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2 Proofs 



Proof of Lemma 1. In [2] we proved that for any exactly-solvable operator T, 
the eigenvalue problem T(p n ) = X n p n can be written as a linear system MX = 
Y, where X is the coefficient vector of the monic nth degree eigenpolynomial 
p n with components a n fi, a„.i, a n ^2, ■ ■ ■ , a, ljn _i, and Y is a vector and M is an 
upper triangular n x n matrix, both with entries expressible in the coefficients 
of T. With T = £j=i Qj-D^', = Y.TJf ] and p„(z) = £™ =0 a n>i z\ the 

eigenvalue A„ is given by 

k . 
\ - ST n - 
2- °W („-,■)!' 

and the diagonal elements of the matrix M are given by 

k 

M l+1 . l+1 = T A " E a « 

l<j<roin(t,fc) 3=1 

for i = 0, 1, . . . , n— 1. The last equality follows since — j)\ = for i < j < k 
by definition (see Lemma 2 in [2] ) . In order to prove that p n is unique we only 
need to check that the determinant of M is nonzero, which implies that M is 
invertible, whence the system MX = Y has a unique solution. Notice that M is 
upper triangular and thus its determinant equals the product of its diagonal el- 
ements. We therefore prove that every diagonal element Mj+i^+i (i 6 [0, n— 1]) 
is nonzero for all sufficiently large integers n for an arbitrary T as above, as well 
as for every n if degQj = j for exactly one j. 



From the expression 



3=1 



(n-j)! (i-j)'. 



it is clear that Mj+i,i+i ^ for every i e [0, n — 1] and every n if otjj ^ for 
exactly one j, that is if degQj = j for precisely one j, and thus we have proved 
the second part of Lemma 1. 



Now assume that deg Qj = j for more than one j and denote by jo the 
largest such j. Then otj j ^ and we have 



3o 



3=1 



(n-io)! 



1 



(n-j)\ 
n\/(n-j )\ 



0- - jo)! 



1<3<30 



i<i<3'o 



- jo) 
n!(i - j 
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The last two sums in the brackets on the right-hand side of the above equality 
tend to zero when n — > oo, since jo > j and i < n — 1. Thus for all sufficiently 
large n we get 



n! 



("-jo)! 



"Jo, Jo I x 



n!/(n-j )! 



^0 



for every i e [0,n — 1], and we have proved the first part of Lemma 1. □ 

To prove Theorem 1 we need the following. If q n is a polynomial of degree 
n we construct the probability measure /j, n by placing a point mass of size ^ 
at each zero of q n . We call /i n the root measure of q n . By definition, for any 
polynomial q n of degree n, the Cauchy transform C n j of the root measure [iffi 
for the jth derivative q^ is defined by 



(n-j)q^(z) J z-C 

for j = 0, 1, ... ,n — 1, and it is well-known that the measure fin ^ can be re- 
constructed from C n j by the formula ^ — ^ • 9C ^' 3 where d/dz — ^(d/dx + 
id/dy). 

Proof of Theorem 1. Let T = X^j=i QjD^ and denote by jo the largest 
j for which degQj = j. Note that since T is degenerate we have jo < k. From 

litii 

P { i\z) 



the definition C n j(z) = — ^ — piy— of the Cauchy transform of Pn\z) we get 

(n-j')Pn (z) 



C n ,o(2)C„ 4 (z) • •■C Tli _ 7 -_i(2) • n(n - 1) • • • (n-j + 1) 



= t — TT c„. m (z). 

n - i ! - L - L 



n! 



With the notation = X^=o ^'.i 2 * we have A„ = X/^Li a j,j (n-j)i ' ^ ow 

dividing the equation T(p n ) — X n p n by p„ we obtain 



Qk(Z) T-r- + yfc-l(^) ^ h . . . + pr = > 



Pn(z) ' ~° K p n (z) ^' K ~'p n {z) j-^ 1,0 {n- j)\ 

, fe-1 , fe-2 

Qk(z) , ' , 11 C n , m (z) + Qfc,!^) — — -— 11 C„, m (z) + . . . 

V m=0 V '' m=0 

Tl x — > Ti' 

• ' • + QlW (t7=T)! C "' oW = Y.^—j jT • ( 2 ) 
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k-l 



Dividing both sides of this equation by ^"'^i we get 

(n-fc)! 1 Q fc -i(z) 

(n-fc + l)!C n , fe _i(z) Q k {z) 

1 Qk-2(z) , 



Qfc(^) JJ Cn,m(z) 
(n-fc)! 



(n - fc + 2)! C n ,fc_i(«)C I , 1 fc_a(«) 



(n-fc)! 



1 



( n - 1 ) ! n^r= 1 iCn,m(2)Qfe(^ 



5> 



(n-fc)! 



(3) 



Now assume that all zeros of all p n are uniformly bounded. Then we can take a 
subsequence {p ni } such that all the corresponding root measures are weakly 
convergent to a compactly supported probability measure. Then all Cauchy 
transforms C ni . m will be uniformly convergent to a non- vanishing function out- 
side some large disc, which in particular contains all the roots of Qk(z). Since 
jo < fc, the right-hand side of tends to zero when n — > oo. On the other 
hand, in the left-hand side of (J3J), all terms in the bracket except for the con- 
stant term 1 tend to zero when n — > oo, and thus the limit of the left-hand side 
equals lim^^oo Qk(z) Y\m=o Cn,m(z) = K 7^ 0, and we obtain a contradiction 
when n — > 00. □ 



In order to prove Theorem 2 we need the following two lemmas, where Lemma 
2 is used to prove Lemma 3. 

Lemma 2. Let z n be a root of p n with the largest modulus r n . Then, for 
any complex number zq such that \zq\ — ro > r n , we have \C n .j(zo)\ > for 
all j ;0. 



f djj?^ (C) _ P { J + 1) (z) 



With £ being some root 



Proof. Recall that C n j(z) := J 

of pn (z) we have \(\ < \z Q \ by Gauss-Lucas theorem. Thus = ^ • i_g/ Zo 
~z~ ' T~H ^ere \6\ = \(/zq\ < 1. With w = — Ly we obtain 



\w - 1| = 
and thus 

\Cn,j(zo)\ 



\0\ 



\9\\w\ < \w\ &\w-l\< \w\ Re(w) > 1/2, 



zq-( 



> 



1 



1 










J 


1-9 


ro 


1 J) (C) 


" 2r J 





1 



□ 



Lemma 3. Let T = V* O,!)' = E*=i (EiS^ ".!■■'') i)j be a degener- 
ate exactly-solvable operator of order fc. Wlog we assume that Qk is monic, i.e. 
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&k,degQ k — 1- Let z n be a root of p n with the largest modulus r n . Then the 
following inequality holds: 



fe-ldcgQ, Jfe-J-degth-H 

i<E E \^ k - j ^ 



E 



\ a k,; 



3=1 t=0 



V ; 0<i<dcgQ fc '« 



dog Q fc -i ' 



(4) 



Proof. From C nJ (z) = ^ )p u*\ z) we get 



P ( n\z) 



V j<k. (5) 



(n — k + l)(n - fc + 2) • • • (n - j) ft^. C n , m {z) 
Inserting z n in the eigenvalue equation Tp n (z) = X n p n (z) we obtain 

fe-l /dogQj , degQfc , 

E ( E "j.^nJPn^n) + ( E a k,i z n)PnH z n) = KPn(Zn) = 0, 
j=l ^ i=0 ' ^ i=0 ' 

and after division by z^ cg< ^ k p n k \z n ) we obtain 
fc-i 



E ( E a J< 1 dcgc 
j=i v i=o z ™ 



dcgQfc-i 7 (fe) 



Pn ^ (Zfi 



E 



a*,- 



i<dcg Qfc 



dog Qfc-i 
Zn 



1 = 0. 



Thus, applying |5} and Lemma 2, we obtain 

I \Vn i \z n ) 



1 



< 



< 



< 



fe-l , degQj 

e(e« 

j=l V i=0 
fe-l degQj 

E E a 

3=1 i=0 
fe-l dog Q, 



dog Q fc -i / (fe) 
A,rt, ' ■ ' 



+ E 



i 



i 



3,1 dog Q fc -i 



Pn (z n ) 0<i<degQ fe 



Q!fc,t- 



E 



dog Q fc -i 



a 



dog Qt-i 



EE 

fc-ldcgQj 

dcgQ 



\Pn \z n )\ 0<i<dogQ fc r 
1 



dog Qfc-i 



(n - k + 1) • • • (n - j) IlmJj- |C„, m (z„) 



E 



i a fc,i 



dog Qfc-i 



0<i<dcgQ, 



E 



3=1 i=0 '' 
fe-l dog Qj 

E E K*i 2 

3=1 »=0 



(n - k + l) k -i — r 

y 1 0<i<dogQfc ' 



\CXk,i\ 
dog Qfc-i 



fe— j n 



P fe-j-dcg Qfc+i 



E 



y ' 0<i<dogQfc 



Qfc,i 
dog Qfc-i ' 



□ 



The proof of Theorem 2 follows from Theorem 1 and Lemma 3. 

Proof of Theorem 2. Applying Theorem 1, we see that the last sum on 



11 



the right-hand side of inequality ^ in Lemma 3 tends to zero when n — > oo. 
Now consider the double sum on the right-hand side of (|4]). If, for given i and 
j, the exponent (k — j — deg Qk + i) of r n is negative or zero, the corresponding 
term tends to zero when n — > oo by Theorem 1. We now consider the remaining 
terms in the double sum, namely those for which the exponent (k—j —deg Qk+i) 
of r n is positive. If r n < co(n — k + l) 7 where Co > and 7 < , _ . , n f° r 
given jg [1,^-1] and given i 6 [0, degQj], then the corresponding term 



(n-k + l) k J V(„ - k + 

in the double sum tends to zero when n — > 00. Thus assume that r„ < Co(n 
k + l) 7 where Co is a positive constant and 7 < 6 where 

+ fc — j + k — j 

b = mm ■ ■ : — — = mm 



jeii k -i] k - j + i - deg Q k je[i,*-x] k-j + deg - deg Q fe ' 

and where the notation min + means that we only take the minimum over posi- 
tive terms (k — j + i — deg Qk) and (k — j + deg Qj — deg Qk) Then every term 
in the double sum tends to zero when n — > 00, and we obtain a contradiction 
to (j4|) when ri — > 00. Thus for all sufficiently large integers n we must have 
r n > co(n — k + l) 7 for all 7 < 6, and hence liminfn^oo > Co for any 7 < 6. 
But for any such 7 we can form 7 = 7 -^ for which 7 < b and 7 < 7 , and thus 
lim n _»oo ^7 = 00 for all 7 < 6. □ 

Proof of Corollary 1. For this class of operators we have 

, + k-j 

:= mm : — — 

36[i,*-i] k-j +degQj - degQ k 

+. k-j k- jo 1 

= mm = = 1, 

j£[i,fc-i] k-j + degQj - j k-jo 

and the proof is complete applying Theorem 2. □ 

Proof of Corollary 2. For this class of operators we have 

k-j 



b := 



+ 
mm 



je[iM-i] \ k - j + dcgQj - degQk 

+ ( k-j \ . k-j k-jo 

mm — = mm — 



j6[i,fc-i] \k - j + degQjJ je[i,j ] k k 



4 On the left-hand side in the expression for b above we take the minimum over 
i 6 [0,degQj], so we can put i = deg Qj in this expression. Thus with b = 
min+ (1 >fc _ 1] k _ J+dc ^Q J j _ dcgQk we get that 7 < k _ j+ k r J dceQk for every j G [1, k - 1] and 
every i £ [0, degQj]. Thenifr n < co(n — + and 7 < 6, every term with positive exponent 
(k — j + i — deg<2 fe ) wm tend to zero when n — > 00. 
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where the third equality follows from choosing any j for which degQ., = j, 
and the minimum is then attained for j = jo (note that for j > jo we get 
( k — J ') j ( k — j ' + deg Q j ) = 1 > (fc — jo ) /k) , and the proof is complete applying 
Theorem 2. □ 

Remark. Note that for the class of operators considered in Corollary 1 the 
Main Conjecture claims that linin^oo ^ = ct for some ct > 0, since d := 
maxjgyo+i^] { jJdJgQ ) — k-jo ~ 1 max ' mum i s attained by choosing any 
j > jo such that degQj = jo, e.g. j = k), and for the class of operators consid- 
ered in Corollary 2 the Main Conjecture claims that lim„^oo ( k -" )/ k — ct for 
some ct > 0, since 

, ( j-jo \ ( 3 - io \ k- j 
a := max : — = max = — ; . 

je[j +i,k] \j -degQjJ je|>o+i,fe] \ 3 J « 

Remark. For a class of operators containing the operators considered in 
Corollaries 1 and 2 we can actually prove that the conjectured upper bound 
linin^oo sup % < c\ implies the conjectured lower bound linin^oo inf % > Co 
where c\ > Co > 0, see Theorem 5 in Section 5.2. 

Proof of Theorem 3. Clearly degQj = jo since there exists at least one 
such j < k. Set 

jo dog Q k 

T = Q jo D*> + Q k D k = J2 a Jo^ l D JO + J2 a k^D k , 

i=0 i=0 

where ctj j ^ 0, and where we wlog assume that Qk is monic. From inequality 
(f4]) in Lemma 3 we have 

J° r i-degQ k +k-j i 

t=0 y ' 0<KdcgQ fc '« 

30 r i-dcg Qk+k-jo 



i=0 



(n-k+ l) fe -io 



where we choose n so large that e < 1 (this is possible since e — > when n — > oo 
due to Theorem 1). Thus for sufficiently large n we get 

Jo „j-deg Q fc +*;-jo 

co < ^2\a j0ti \2 k -^ " 



(n-k + l) fc -io 
i=0 v ' 

Jo fe-degQ fc 

< > Irv - \2 k ~ jo — 

^ 2J a *>,*M („_fc + l)k-io 

r fc-dcg Q fe 

= if- 



(n - k+ l) fc -J'o ' 
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where 1 — e = cq — » 1 when n — > oo, and if > since (Xj j ^ (the second 
inequality follows since i < jo). Thus 

for sufficiently large integers n, and hence there exists a positive constant c = 
(jyjQi/CA-degQ*) such that 

lira inf — ? r > c. 

n— >oo ( k -io ) 

fi \ fc-dcg Qfc / 

Finally, it is clear that for this two-term operator 

3- jo \ k - 3o 



d 



je[? +i,fc] \j - deg Qj ) k- deg Q k 
and we are done. □ 

Remark. If, in Theorem 3, Qk is a monomial (i.e. Qk — z dcs ® h ), then 
there exists a positive constant c such that r n > c(n — k + l) d for every n, 
where d := max je{jo+lik] (j^q-) = k - A J° Qk ■ Th is is easily seen from the 
calculations in the proof of Theorem 3 since 2o<i<degQjs \ a k,i\ A» S Q h -* on the 

k — deg Q k 

right-hand side of Q vanishes, and therefore 1 < K ^-k+i) k -io for every n. 
From the second part of Lemma 1 we know that for this class of operators there 
exists a unique eigenpolynomial p n for every n, and the conclusion follows. 

Proof of Theorem 4. For this class of operators (j — degQj) > (k — degQk) 
for every j > jo and thus 

j -jo \ It- 3o 



d := max . 

ie[io+i,fc] \j - degQj J fc - degQ* 

Assuming that Q^ is monic we have the inequality 

i<y y \ aji \2 k -*i±— — - k —+ y } a ™ . (6) 

3 = 1 t=0 v ; 0<i<dcgQ fc '« 

by Lemma 3. The last sum here tends to zero when n — > oo by Theorem 1. 
Considering the double sum on the right-hand side of ((6]) we see that for every 
j we have (since i < degQj) that 

degQj k _ j+i _ degQk degQj fc-j+degQj-degQ* 

y | Q |2*-jI»_ . — __= y | aji |2 fc -^ — - — r,r dcgQj 



n fc-j+degQ 3 — degQn / \ 

— (2^|a J - degQ .|+ y 2* '<-,, 'V desQ J 

i<dcgQ., 



(n - k + l) k -i 

fc-j+deg Qj-degQfe 



~ Kj ' n n {n-k+l) k -i ' (7) 
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where 



K j>n = 2 k J \a jAcgQj \+ 2 * " 1° 



«<dcg Qj 

where ifj,rj > since a^degQ* 7^ 0. Also, K^ n < 00 when n — *■ cxd since 
(i — dcg Qj) < for the exponent of r n , and then using Theorem 1 (note that 
Kj :n — > 2 fe ~- J loij^cgQj I when n — > 00). With the decomposition 



A = {] 
B = {j 
C={j 



degQj =j}, 

deg Qj < j and (fc - j + deg Qj - deg Q k ) > 0}, 
deg Qj < j and (fc - j + deg Qj - deg Q k ) < 0}, 



and using J7J) we see that inequality © is equivalent to: 



v^V"^ fc- ■ r k l ^ ]+i ^ dcgQk \ - lafc. 



3=1 i=0 v ' 0<i<degQ k ™ 

„Ai-dcgQ fe k-j+dcgQj-degQk 

2^ ( n -k+ l) k -i ^ 3n (n-k + l) fe -J 

fc-j+dogQj-degQfc i i 



(n-fc + l) fe -^ ^ dcgQ fc - 4 - 

jeC V ; 0<i<degQ k 

The last two sums on the right-hand side of this inequality both tend to 
zero when n — > 00, the last one due to Theorem 1, and the sum over C since 
(j - degQj) >{k- degQfc) <S=> (fc - j + degQj - degQ k ) < for every j e C 
by assumption, and then applying Theorem 1. Therefore, when n —> 00, we get 
the inequality 

r fc-dcgQ fc r k-j+degQj-degQ k 



where 



'(n-fc + l) fc -J ^ J '™ (n-fc + l) fc -J 
jeA ' jeB y 1 



fc- j +dcg Q j -dcg Q k , 
_ , \ -» „ \ "* \Otk,i 



(n-fc+l) fc -J ^ dc g Q fc - l - 

Note that cq — > 1 when n — > 00. 

Now assume that B is empty. This corresponds to an operator such that 
(j — degQj) > (fc — degQk) for every j for which degQj < j. Then inequality 
© above becomes 



co < XX>7— n 



(n-k + l) k ~i 

r fc-degQ fc / l 

= (n- n fc+l)*-Jo ^'•"(n-fc + l) J o- J 

„fc-deg Q k 
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where K~a is a positive constant which is finite when n — > oo, since jo — j > 
for every j E A\{jo} (recall that jo is the largest element in A by definition). 
Thus for all sufficiently large integers n we have 

r n >(^) 1/{k - dCsQk \n-k + l)-^, 

ft A 

and therefore there exists a positive constant c = (1/ Ka) 1 ' ^ k ~ des< ^ h ' such that 

lim inf — ; r- > c, 

n->oo ( k -JQ ) 

n \k-d B g Q k I 

and we are done. 

Now assume that B is nonempty. Again inequality (jSJ holds, i.e. 

r fc-dogQ fc jk-i+degQj-degQk 
CO < J2 K i* ( n " k + + J2 K >> "(„ _ k + i 

where c — ► 1 when n — > oo. From ([5]) we have 

E fc-dcgQfc fc-dcgQ fc 
If. " < If,, " 

(n - k + \) k -i ~ A (n-k+ l) k -i° 

jeA v ' x ' 

for the sum over A for large n and thus 



Co 



r k~dagQ k k-j+degQj-desQh 

S ^ 3 ' n (n - k + l) k -> ^ J ' n (n-k+l) k -i 

„fc-dcgQ fc fc-i+dogQ 3 -degQ fc 

" A (n - k + l) k -i« + f^ B °' n (n-k+ l) k -J 



fe-dogQfc / r ^gQj-j \ 

jeB 



(n-k + l) fc -io 

r k— deg Q k 



< K AB 



(n- k + l) k -i° ' 



where Kab is a positive and finite constant when n — ► oo (note that Kab ~> Ka 
when n — > oo, since (deg Qj — j) < and j — j > for every j S £>). Thus for 
all sufficiently large integers n we have 

r„>(-^) 1/(fc - dcsQfc) (n-fc + l)^:, 

ft-AB 

so there exists a positive constant c = (l/J ; GiB) 1/ '( fc ~ dos( 2' 1 ) such that 

lim inf — -, r- > c. 

n— >oo ( k -JQ ) 

n\ fc — deg Q k J 

□ 
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3 Main Conjecture and its Corollary 



In this section we explain how we arrived at Main Conjecture (see Section 1) 
and obtain as a corollary of our method the (conjectural) algebraic equation 
satisfied by the Cauchy transform of the asymptotic root measure of the prop- 
erly scaled eigenpolynomials. 

How did we arrive at Main Conjecture? 

Let T = QjDi = J2j=i ( 2i=o a j,i z% )^ >: ' ^ e an arbitrary degenerate 

exactly-solvable operator of order k and denote by jo the largest j for which 
degQj = j. Wlog we assume that Qj is monic, i.e. otj Qi j a = 1. Consider the 
scaled eigenpolynomial q n (z) = p n (n z), where p n (z) is the unique and monic 
nth degree eigenpolynomial of T, and d is some real number. The goal is now 
to obtain a well-defined algebraic equation for the Cauchy transform of the 
root measure fj, n of the scaled eigenpolynomial q n when n — ► oo, and as we will 
see in the process of doing this, we are forced to choose d as in Main Conjecture^ 

Basic assumption. When performing our calculations we assume that the 
root measures , M» ^ > A*n^ ■ ■ ■ , fMi ^ of the scaled eigenpolynomial q n (z) and 
its derivatives up to the fcth order exist when n — > oo and that they are are 
all weakly convergent to the same asymptotic root measure Thus the cor- 
responding Cauchy transforms are all asymptotically identical, and we define 
C(z) :— lim„_ >oo C n j(z) for all j G [0, k—1], where C(z) is the Cauchy transform 
of fj, and is considered for z's away from the support of [i. Computer exper- 
iments strongly indicate that this assumption is true - for details see Section 4.2. 

From the definition of the Cauchy transform we obtain 

n^w = n ; qn , ; i 

4 =o {n-j)qn[z) 



i=0 



q^{z) $\z) q£\z) 



nq n (z) (n-l)qW(z) ( n - 2'jq^ (z) 
qir 1] (z) $\z) 



(n-j + 2)q^- 2) (z) (n-j+ l)q { ri- 1} (z) 



n(n - 1) • • • (n - j + l)q n {z) ' 



5 It is already well-known that for the Laguerre polynomials, which appear as eigenpolyno- 
mials for a second order exactly-solvable operator, the largest root grows as n when n — > oo 
and thus d = 1 in this case, which is consistent with Main Conjecture. 

6 Conjecturally supp /i is a tree, see Section 6 on Open Problems. 
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and thus our basic assumption implies 

C 3 (z) = lim Y[C nii (z) = lim g " 1 ; (10) 

n— >oo n-»oo n(n — 1) • • • (n — j + l)q n {z) 

In the above notation consider the eigenvalue equation Tp n (z) = \ n p n {z), 
where the eigenvalue A„ is given by 

k | jo j jo 

A « = E "to r „"',- v = E a « f—fu = E a w n ( n - 1) • • • (n - i + 1). 
3=1 v 3=1 v J> ' j=i 

Clearly this sum ends at jo since otj j = for all j > jo by definition of jo as 
the largest j for which deg Qj — j. We then have 

Tp n (z) = \ n p n (z) 

<^ 

fe .dcgQj j 

E ( E a 3,i Z )Pn ) (z) = E a 3d n ( n -!)•••(«- j + l)Pn(z). 
j=l \ l= ' 3=1 

Substituting z = n d z in this equation we obtain 

E ( E )i#V*) = E - 1) • • • (n - j + l) P „(« d Z ), 

3=1 »=o ' 3=1 

and with q n (z) = p n {n d z) we get 

E( E ^ ndU-i) )g» = E ~ : ) ' ' ' ( n ~ J + 

3=1 «=0 ' 3=1 

Dividing this equation by ^"-. o)! g n (z) = ra(n - 1) • • • (n - j + l)q„(z) we get 



k / dc sQj N „(j) 

E E a ^ 

j=l v i=0 



'* n d ^ ) n(n - 1) • • • (n - j + !)«»(«) 



A n(n~ l)---(n- j + 1) 
" 2-""„„ ( L)...( n j • if 

Consider the right-hand side of (fTTj) . Since j < j all terms for which j < j (if 
not already zero, which is the case if otjj = 0, i.e. if degQj < j) tend to zero 
when n — > oo, and therefore the limit of the right-hand side of (fTTj) equals 

A n(n - 1) • • • (n - j + 1) 

lim > a, 7 — ; : ; = Q.j n >„ = 1, 

n-ooj-^ 3 ' 3 n{n-l)---(n- jo + l) 30,30 
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since we assumed that Qj is monic. Now consider the jth term in the sum on 
the left-hand side of (fTTjl . It equals 



deii 



i=0 



TV 



n(n-l)---(n-j + l)q n (z) 



de j% 3 a _j_ q { n 3) (z) n(n- !)■■>- j + 1) 

^ n d 0'-*) n(n - 1) • • • (n - j + l)g„0) n(n - 1) • • • (n - jo + 1) 

ST „ z Tin I \ n{n-l)---{n-j + l) 



^^-i) »(n-l)...(»-jo + l) 

_ Y[ r ( \ n(n-l)---(n-j + 1) n>° 

- 2^ a 0,i d{j - i)+jo -j • 11 ^n,iW ^ „( n _l)...( n _j +!)• 

Taking the limit and using the basic assumption (|10j) we obtain 

a, 



lim \ 

n — >oo ^ — ' 



__ n '^n^'-') n(n-l)---(n-j + l)q n (z) 
lim > a,-j— jtt — — tCHz) 

i=0 



for the jth term and thus, taking the limit of the left-hand side of (fTTj) we get 

fc / dCg9j 7 l \ n {j) (r\ 

iim y f y a z » g - (z) 



v , ,, ' n d O--i)y n (7i-l)---(n-jo + l)9„(2;) 
= limyf y aji ... *. . )CUz). 

j=l v i=0 7 

Adding up, the following equation is satisfied by C(z) for z's away from the 
support of fi: 

j=l v i=0 7 

In order to make (|12[) a well-defined algebraic equation, i.e. to avoid infinities 
in the denominator when n — > oo, we must impose the following condition on 
the real number d in the exponent of n, namely 

d(j-i)+jo-3>0 d>^— ^ 
for all j £ [l,fel and alii G [0,degQ,l. Therefore we take d = max ssim t 2 ^-), 

i6[0,degQ 3 -] 

but this maximum is clearly obtained for the maximal value of i for any given 
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j, so we may as well put i = degQj. Our condition then becomes d = 
maXj- e [i J k] ( jJfcgQ ) ' anc ^ clearly the maximum is taken only over j for which 
Qj(z) is not identically zero. Finally we observe that since T is degenerate we 
have jo < k and thus we need only take this maximum over j 6 [jo + l>k], 
since there always exists a positive value on d for any operator of the type we 
consider. Thus our condition becomes: 

, ( j — jo \ 
d = max — — . 



Corollary of Main Conjecture. 

In the above notation (recall that ay j = 1 by monicity of Qj ) , the following 
well-defined algebraic equation follows immediately from inserting d as defined 
in Main Conjecture into equation (|12[) and letting n — > oo: 

Corollary. The Cauchy transform C(z) of the asymptotic root measure fj, of the 
scaled eigenpolynomial q n (z) — p n (n d z) of an arbitrary exactly- solvable operator 
T as above satisfies the following algebraic equation for almost all complex z in 
the usual Lebesgue measure on C: 

z*>C*>(z) + a lAcgQ] z d ^C'3(z) = 1, 
jeA 

where A is the set consisting of all j for which the maximum 

d := maxj e y 0+1)fc ] ( jJ^° Q . ) is attained, i.e. A = { j : (J - jo) /{j -degQj) = d}. 

For detailed calculations see Section 5.1. 

4 Numerical evidence 

4.1 Evidence for Main Conjecture 

In the table on the last page of this section we present numerical evidence on 
the growth of r n — max{|z| : p n (z) = 0} which supports the choice of d in Main 
Conjecture. We have performed similar computer experiments for a large num- 
ber of other degenerate exactly-solvable operators, and the results are in all cases 
consistent with Main Conjecture. Next we present some typical pictures on the 
zero distribution of the appropriately scaled eigenpolynomial q n (z) — p n (n d z) 
for some degenerate exactly-solvable operators, where p n denotes the nth degree 
unique and monic polynomial cigenfunction of the given operator T. Conjec- 
turally the zeros of q n {z) are contained in a compact set when n — » oo. 
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Fig.l: T 1 = zD + zD 2 + zD 3 + zD 4 + zD 5 . 




roots of roots of roots of 

Q5a(z) =_p 5 o(50z) q 75 (z) = p 75 (75z) qwo(z) = Pioo(lOOz) 

Fig. 2: T 2 = z 2 D 2 + D 7 . 




roots of roots of roots of 

q 50 (z) = P 5 o(50 5 / 7 z) q 75 (z) = P75 (75 5 / 7 z) q 100 (z) = Pwo {l00 5 / 7 z) 

Fig. 3: T 3 = z 3 D 3 + z 2 D A + zD 5 . 



roots of roots of roots of 

«5o(*) = P5o(50 1/2 z) q 75 (z) = pn>V& ,2 z) qioo(z) = proo^O 1 / 2 *) 

4.2 On the basic assumption 

Let us show some examples of pictures supporting the basic assumption, namely 
that the root measures of q n (z) and its derivatives up to the fcth order exist when 
n — > oo and are all weakly convergent to the same measure \i. 

Fig. 1: T 7 = zD + D 3 and q n (z) = p ri (n 2 ' 3 z). 




roots of qioo(z) roots of q' wo (z) roots of qi 00 (z) roots of q'lboi 2 ) 
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Fig. 2: T 9 = zD + zD 4 + z 3 D 7 and q n (z) = p n {n 3 / 2 z) 





roots of qioo(z) roots of q'ioo( z ) roots of q" ao (z) roots of qio (z) 





roots of q[™Q (z) roots of q^ (z) roots of q^ (z) roots of q^Q (z) 



4.3 On the Corollary of Main Conjecture 

The algebraic equation in the corollary of Main Conjecture satisfied by the 
Cauchy transform of the asymptotic root measure of the scaled eigenpolynomial 
indicates that the asymptotic zero distribution depends only on the term z JO D JO 

is at- 



j-dcgQj 



and the term(s) Oj,dcgQj z dcs< ^ j Z3 J for which d = max 3e [ 3o+li t] ( 

tained0 Thus any term QjD^ in T for which j < jo or such that (j — ja)/{j 
degQj) < d is (conjecturally) irrelevant for the zero distribution when n — > oo. 



To illustrate this fact we now present some pictures of the zero distributions 
of the scaled eigenpolynomials for some distinct operators for which the Cauchy 
transforms C of the corresponding scaled eigenpolynomials q n satisfy the same 
equation when n — > oo, namely the equation in the Corollary of Main Conjec- 
ture in Section 3. 



As a first example consider the operator T4 



z 3 D 3 + z 2 D 5 . 



Here d 



max jebo+i,fc] 



3 -Jo \ — 



(5- 
Pn{n 2/3 



3)/(5 — 2) = 2/3, the corresponding scaled eigen- 



and we have z 3 C 3 + z 2 C 5 — 1 for the Cauchy 



^J-degC 

polynomial is q n (z) - 
transform of q n when n — > 00. Now consider the slightly modified operator 
T 4 = z 2 D 2 + z 3 D 3 + zD A + z 2 D 5 + D 6 and note that d is obtained again (only) 
for j = 5 (for j — 4 we have (4 — 3)/ (4 — 1) = 1/3 < 2/3 and for j — 6 we 
have (6 — 3)/(6 — 0) = 3/6 = 1/2 < 2/3). We therefore obtain the same Cauchy 
transform equation as for T4, and hence the terms z 2 D 2 , zD 4 and D 6 in T4 can 
be considered as irrelevant for the zero distribution for sufficiently large n. The 
pictures below clearly illustrate this. 



7 Recall that we normalized T by letting Qj be monic, i.e. ctj Ol j = 1. 
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T 4 = z 3 D 3 + z 2 D a , T 4 = z 2 D 2 + z 3 D 3 + zD 4 + z 2 D 5 + D e , 

roots of qioo(z) = pioo(100 2/,3 z) roots of qwo(z) = pi o(100 2/3 z) 



However, instead of D 6 , we may add the more "disturbing" term zD 6 to T4. 
Note that for the operator T 4 = z 2 D 2 + z 3 D 3 + zD 4 + z 2 D 5 + zD 6 for j = 6 
we have (6 — 3)/(6 — 1) = 3/5 = 0.6 < 2/3. Adding any term QjD^ such that 
U " jo)/ (j — deg Qj) < d to a given operator, it is clear that the closer the value 
01 (.? ~ Jo) / (.7 — degQj) is to d (in this case 2/3), the more disturbing it is in 
the sense that it requires larger n for the corresponding zero distributions to 
coincide. See pictures below. 



T4, roots of T4, roots of T4, roots of 

giooC*) = Pioo(100 2 / 3 z) q wo (z) = Pl0 o(W0 2 / 3 z) q 60 o(z) = p 600 (600 2 / 3 z) 

Increasing n however, experiments indicate that the zero distributions of the 
scaled eigenpolynomials of T4 and T4 coincide, as they (conjecturally) should. 

As a second example, consider the operators T5 = z 5 D 5 + z 4 D 6 + z 2 D 8 and 
T 5 = z 2 D 2 + z 5 D 5 + z 4 D 6 + zD 7 + z 2 D 8 , whose scaled eigenpolynomials q n (z) = 
Vn(yi}l 2 z) both satisfy the Cauchy transform equation z 5 C 5 + z 4 C 6 + z 2 C 8 = 1 
when n —> 00. In the pictures below we see that the terms z 2 D 2 and zD 7 of T5 
seem to have no effect on the zero distribution for large n. 



T 5 = z 5 D 5 + z 4 D 6 + z 2 D 8 , T 5 = z 2 D 2 + z 5 D 5 + z 4 D e + zD 7 + z 2 D 8 , 
roots of qioa(z) = pioo(100 1/2 z) roots of qwo{z) = pioo(100 1//2 z) 
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Finally note that for j and for any j for which d is attained, it is only 
the highest degree term etj,degQj z dce ®i of Qj that is involved in the Cauchy 
transform equation. Consider for example the following case, where adding 
lower degree terms to aj,degQj z dog< 3 J ' in the (relevant) Qj seems to have no 
effect on the zero distribution for large n. Below T§ = z 3 D 3 + z 2 D 6 , and 
T 6 = [(1 + 13i) + (24i - 3)z + Uiz 2 + z 3 ]D 3 + [(22* - 13) + (-9 - Ui)z + z 2 }D 6 . 
Note the difference in scaling between the pictures. 



X 




T$, roots of 
9ioo (*) =Pioo(100 3 / 4 z) 



T 6 , 

qwo(z) = 



roots of 
p 100 (100 3 / 4 z) 



T 6 , roots of 
gsoo (z) =p 5 oo(500 3 / 4 z) 
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Operator 


n 


r n experimental 


r n conjectured 


Ti = zD + zD 2 + zD 3 + zD 4 + zD 5 


50 


„ _ rr ,(] Utf7>ig& 

2.7 • 50 u aD ' oao 


r „i 
ci • 50 1 


100 


2.7 • 100 


ci • 100 1 


200 


o -7 onnU 992657 

2.7 • 200 


ci • 200 1 


250 


2.7 • 250 


Ci • zol) 


T 2 = z 2 D 2 + D 7 


50 


1.3 • 50 U b7iy77 


1-^^/7 

c 2 • 50 5/7 


1UU 


1 1 1 nnO. 694847 
l.o • 1UU 


C2 * 1UU ' 


200 


1 Q onnO. 706226 
1.3 • 200 


onn5/7 
C2 • zUU ' 


onn 
300 


1 i qnnO. 710085 
1.3 • 300 


C2 * ol)l) ' 


/inn 
4UU 


1 Q /innO. 712043 
1.0 • 4UU 


^ /i nn5 / 7 
C2 * 4UU ' 


T 3 = z 3 D 3 + z 2 D 4 + zD b 


50 


4/3 • 50 u - 4byou7 


c 3 • 50 1/2 


1 nn 
1UU 


A /I 1 nnO. 484824 

4/o • 1UU 


C3 • 1UU ' 
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200 


/I /q onnU. 492832 
4/o • 200 


C3 • 200 ' 
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OUU 


a/i inn - 495592 

4/ o • OUU 


inn 1 / 2 

C3 • oUU ' 


Ann 

4UU 


44/ O * 4UU 


„ . /inn 1 / 2 

03 • ^UU ' 


T A = z 3 D 3 + z 2 D 5 


50 


1 4 5Q0-633226 


cn2/3 

C4 • 50 10 


i nn 

1UU 


i 4 innO. 652141 
±.4t • 1UU 


„ inn 2 / 3 

C4 • 1UU 1 


9nn 

ZUU 


1.4 • zuu 


„, 9f|f|2/3 
C4 • ZUU ' 


inn 

OUU 


1.4 • oUU 


r- inn 2 / 3 

C4 * r3UU ' 


Ann 

4:UU 


1.4 • 4UU 


„, /inn 2 / 3 

04 • ^UU ' 


T A = z 2 D 2 + z 3 D 3 + zD 4 + z 2 D 5 + D 6 


50 


1 a cnO 632811 

1.4 • 50 U DOZO±1 
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1 a inn - 651960 

1.4 • 1UU 
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C4 • 1UU ' 


200 


1 A onnO. 661332 
1.4 • zUU 


C4 • 200 ' 
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OUU 


1.4 • OUU 


r, . inn 2 / 3 

C4 * OUU ' 


/inn 

4UU 


1.4 • 4UU 


n, /inn 2 / 3 

C4 • 4UU ' 


T 5 = z 5 D 5 + z 4 D e + z 2 D s 


50 


i r rnf) 462995 

1.5 • 50 u -^ yyo 


C 5 • 50 ' 


i nn 

1UU 


1.0 • 1UU 


„_ inn 1 / 2 

C5 • 1UU ' 


9nn 

ZUU 


1.0 • zuu 


-, K . onn 1 / 2 

C5 • ZUU ' 


inn 

OUU 


1 ^ inn"- 494 *! 4 

1.0 • OUU 


„ K . inn 1 / 2 

05 • OUU ' 


/inn 

4UU 


1.0 • 4UU 


„ /inn 1 / 2 

C5 • 4UU ' 


T 5 = z 2 D 2 + z 5 D 5 + z 4 D e + zD 7 + z 2 D s 


50 


1.5 • 50 { > /±OMV1 


c 5 • 50 ' 


i nn 

1UU 


1.0 • 1UU 


^ 1 nn 1 / 2 

C5 • 1UU ' 


onn 
zuu 


1 ^ OHO - 491129 
1.0 • ZUU 


?r onn 1 / 2 

C5 • ZUU ' 


inn 

OUU 


1 ^ inn" ■ 494342 

1.0 • OUU 


%. inn 1 / 2 

C5 • r3UU ' 


/inn 

4UU 


1.0 • 4UU 


%. /inn 1 / 2 

C5 • 4UU ' 


T 6 = z 3 D 3 + z 2 D 6 


50 


1 A r nO 7021 1 7 


rr»3 /4 

c 6 • 50 3/4 


100 


1.4 • l00 u - Y2&Yi& 


c 6 • 100 a / 4 


200 


1.4 • 200 u - 73754i 


c 6 • 200 a / 4 


300 


1.4 • 300 U Y4iei4 


c 6 • 300 a / 4 


400 


1.4 • 10U : ' ' 


c 6 • 400 a/4 


T 6 = [(1 + 13i) + (24i - 3)z + \\iz 2 + z 3 ]D 3 
+[{22i - 13) - (9 + Ui)z + z 2 ]D e 

25 


50 


1 4 • 50°- 769260 


c 6 • 50 a/4 


100 


1 4 • lOO 760399 


c 6 • 100 3 / 4 


200 


1.4 • 200 u - 75t)lt > 1 


c 6 • 200 a / 4 


300 


1.4 • 300 u - 7545yo 


c 6 • 300 a / 4 


400 


1.4 • 400 U 7537b5 


c 6 • 400 a ^ 4 



5 Appendix 

5.1 Arriving at the Corollary of Main Conjecture 

The algebraic equation in the corollary of Main Conjecture follows immediately 
from inserting d as defined in Main Conjecture into equation (|12p in Section 3 
and letting n — > oo. If we put d into (|12|) we namely get 

k /dcgQj t . 

E E °w 7 z — v- H z) = '• (13) 

Denote by iV}^ the exponent of n in (| 13[) for given j and i. Thus 
N jti = max _ ( . 3 J (j jo - j. 



j£b'o+i,fc] VJ ~ degQj 

The terms in (|13p for which this exponent is positive tend to zero as n — ► oo. 

First we consider j for which degQj = j, and denote, as usual, by jo the 
largest such j. If j = jo, then i < deg Q j0 = jo and thus for j = jo and i = jo 
we get 

j - Jo 
- degQ 

j - jo 



N j ,h = max , , - — i — 7T (3 ~ *) + io - j 
je|>o+i,fc] V./ - dcgQj 



ie[jo+i,fc] \j - degQ, 
and for j = jo and i < jo we have 



Oo - jo) +jo - jo = 0, 



N 3o,i = max , , ( ■ J 1 ) U + JO - 3 

j - Jo 



> max , ~ — a — FT (-7° ~ -?°) + 3o ~ 3o = 0- 

Thus Nj j n = and A^- 0i i > for i < jo, and for the term corresponding to 
j = jo in (pp we get 

■Jo j 

E a io.' 7 i- J0 ^ ; Ci °^ a jotjo zJ°Ci°(z) = z^C^(z) 

when n — > oo and using Oj- O j = 1. 

Now let j be such that degQj = j and j < jo- Then i < deg Qj = j and 
N~ j = max ( . - ) (j — i) + jo - j 

= P ax , ( ■ 3 a ""In J (■?' ~ •?') + 3o - j = jo - j > 0, 
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and for i < j we get 



J - Jo 



Nj i = max - — - — — (j jo - j 

je[j +i,k] \j - degQj/ 



e[io+i,fc] \j - degQj 



that is Nj : i > for all j < jo such that deg Qj — j and for all i < j. Thus for 
the corresponding terms in (fT3|) we get 

dcg Qj i 

£ E °« ? 1 30 v -u- CJ(z) 



when n — > oo for every j < jo for which deg Q j = j 



3 -Jo 



Now denote by j m the j for which the maximum d = max J(E [ JO +i,A;] ( jJdegQ- 
is attained. Note that there may be several distinct j for which this maximum 
is attained @ Then 

N 3m,degQ jm = . max ( — ^— ) (j - i) + jo - j 



jebo+i.fc] \j - degQj 

Jm - jo 



j m - deg Qj 

jm - jo + jo ~ jm = 0, 



{j m - deg Q Jm ) + jo - jr. 



and for i < deg Qj m we get 

N im,i = P ax , ( ■ J A 3 °^ J (j - + Jo - j 

> / Jm J ° j ( j m - deg Q Jm ) + j - jm 
\j m - deg Q jm J 

= jm - jo + jo ~ jm = 0, 

i.e. Nj mt< \cgQ jm = and Nj mi i > for i < deg Qj m , and for the term corre- 
sponding to j = j m in (fl3l) we get 

E" °w 1 & r, — ; C3m{z) -> a ^Q^ zdc * QomC3 " 



8 Consider for example the Laplace type operator (that is with all polynomial coefficients 
Qj linear) T = zD + zD 2 + . . . zD k . Here jo = 1 and the equation satisfied by the Cauchy 
transform of the asymptotic root measure of the scaled eigenpolynomial q n {z) = p n (nz) is 
given by zC(z) + zC 2 (z) + . . . zC k (z) = 1, since d = maxj 6 [ 2 ,fe] ( j -JdegQ ) = 1 i s attained 
for every j = 2, 3, ... k. 
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when n — ► oo. In case of several j for which d is attained, we put A = {j : 
(j — jo)/(j — degQj) = d}, and for the corresponding terms in lfl3)) we get 

deg Qj j 

E E °W /*-,„ )r — : CHz)^^2a j , dmQj z d ^C^z) 

jeA i=o „ max jE|jo+i,M ( j-dcgQj JO-»)+Jo-j ieA 

when n — > oo. Now consider the remaining terms in (I13p . namely terms for 
which j < jo such that degQj < j, terms for which jo < j < j m , and terms for 
which j m < j < k (clearly this last case does not exist if j m = k) . 

We start with j < jo such that degQj < j. Then i < deg Qj < j and 

N 3,i = P ax , ( ■ 3 - ""V, J U ~ i) + Jo - j 

> max , ( ■ J a J °n ) C? ~ j) + 3o -3= Jo ~ j > 0, 
jebo+i,*] \J ~ degQjJ 

and for the corresponding terms in (| 13j) we have 

deg Qj i 

E E «w r "_ J0 v — - o 

j6{j<Jo:degQj<j} i=0 n lJ ° T cteg / 

when n — > oo. 

Now assume that j m < fc and consider j m < j < k. Clearly j m > jo since 
the maximum is taken over j e [jo + l,fc], and therefore i < degQj < j for 
jm < j <k. Also, 

, J - Jo \ / Jm - Jo \ ^ / J - Jo 
max — = — > 



je[j +i,fc] \ j - deg Q 3 J \j m - deg Q Jm J \ j - deg Qj J ' 
since the maximum is attained for j m by assumption. Thus we get 

J ~ JO \ / . -\ . ■ • I jm — JO 



N jt i = max - — - — — (j - i) + jo - j = — (j - i) + Jo - j 

jeUo+iM \j - deg Q,- / \j m - deg Q jm J 

> ( ■ 3 ~A k n )(j-i)+jo-j>( . 3 ~ j ° n )U-degQ j )+j -j 
\j-degQjJ \j-degQjJ 

= j - jo + jo - j = 0, 

i.e. Njj > for every j m < j < k and every i < degQj. For the corresponding 
terms in (fl3|) we therefore get 

dog Qj i 

E E <*.<— — t4^ ^ — . ^)- 

jm<j<k 1=0 n 



,nax j e [ io + 1 • fc l ( j-dogQj ) C? -*) +Jo - J 



asn-* oo. 
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Finally we consider j < j < j m . Note that this also covers the case j mi < 
j < 3m 2 where the maximum d is attained for both j mi and j m2 . Since i < 
degQj < j we get 

N j4 = max I , 3 30 I (j -i)+j -j= ( , 3m 30 I (j - i) + j - j 
jebo+i,fe] \J - deg Qj J \j m - deg Q jm J 

> ii^k) - j) + J " - J £ ( j^fe) - deg *> + J ° - 3 

= j - jo + jo - j = 0, 

i.e. JVj < > for every jo < j < jm and every i < degQj- Thus for the 
corresponding terms in (jT3j) we get 

deg Qj j 

E E v ; CJ » - 

when n — > oo. 

Adding up these results we finally get the following equation by letting n — > cxd 
in equation |T3|): 

**><7*>(z) + £ a,,dc g Q^ dcg% ^(z) = 1, 

where jo is the largest j for which degQj = j, and A is the set consisting 
of all j for which the maximum d — max_ 7 ' E y _|_i i fc] ( jz^^q~. ) is attained, i.e. 
A={j:(j-jo)/(j-degQ j ) = d}. 

5.2 Theorem 5 

Here we prove that for a class of operators containing the operators considered 
in Corollaries 1 and 2 the conjectured upper bound linin^oo sup(r n /n d ) < c\ 
implies the conjectured lower bound lim JWOO inf(r„/n d ) > cq for some constants 
ci > Co > and where d is as in Main Conjecture. This follows automatically 
from inequality (j4|) in Lemma 3. We have the following 

Theorem 5. Let T be a degenerate exactly- solvable operator of order k which 
satisfies the condition 

. + f k-j \ ( j-jo \ 

b := mm = max =: a, 

ie[i,fe-i] \k - j + degQj - degQ k J je\j +i,k] \j-degQjJ 

where the notation min + means that the minimum is taken only over positive 
values of (k—j+degQj— degQk)- Assume that the inequality r n < ci(n— k+l) d 
holds for some positive constant C\ for all sufficiently large n. Then there exists 
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a positive constant cq < c\ such that r n > cq(u — k + l) d holds for all sufficiently 
large n. Thus 

lim sup — ^ < ci => lira inf — -? > Co- 
Proof. From inequality ([¥]) in Lemma 3 we have 



fc-l deg Qj 

i < e e 1^1^- ^:^^ + E 



r fe_ j+i _ degQfc |q M | 
(n-)fc + l) fc -^ ^ r d °s Qk-i 



.7=1 i=0 v ' 0<4<dogQ fc '« 

fc-l fc-j+dogQj-dcgQfc I I 

< yx/-^ — - — ^+ y ^4 (w) 

] = 1 y ^ ' 0<i<i fc '« 

where the Kj are positive constants. The second sum on the right-hand side of 
(|14p tends to zero when 77 — > 00 due to Theorem 1. We now decompose the first 
sum on the right-hand side of (fl"4|) into three parts. Namely, let 

A = {j : k _ J+dc gQ j _ dcgQk = d}, and note that (k — j + degQj - degQ k ) > 
here since d > 0. 

B = {j : k _ J+dc gQ 3 ._ dcgQk > d}, and note that (k - j + deg Qj - degQ fe ) > 
here since d > 0. 



C = {j : (k -j + degQj - degQfe) < 0}, and note that j < k in l[T4")l. 

Clearly due to the condition b = d there are no terms for which -, — . , , fc Z J , — ^- < 

J fc-j+dcgQj-degQfc 

d and (k — j + deg — deg Q k ) > both hold . 
If j € A then 

jJc-j+dcgQfc-dcgQfc / r \ fe-i+dcgQj-dogQfc 



(77 - + l) fe -J V(77-fc + l) d 

for the corresponding terms in the sum on the right-hand side of |T 



If j £ B then d(k — j + degQj — degQ k ) < (k — j), and this inequality to- 
gether with the upper bound r n < c x (n — k + \ ) d which we assume holds for all 
sufficiently large 77, gives us 

r fe-i+degQ 3 -degQ fc Ci (n - fc + \)d(k- j+dcgQj -deg Q k ) 



(77- fc+ 1)^' ~ (77-fc+l) fc ^' 

when 77 — > 00 for the corresponding terms in (Q2 



If j S C then (k — j + deg Qj — deg Q k ) < and we get 

fc- j+dcg Qj -deg Q k 
— ► 

(77 - k + l) k ~-> 
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when n — ► oo for the corresponding terms in (TT4")) due to Theorem 1. Note that 
if (k — j + deg Qj — deg Qk) — the corresponding term tends to zero when 
n — > oo since j < k in (|14p . 

With this decomposition of the first sum on the right-hand side of the last 
inequality in (114|) we can write 

fe-l k-j+dcgQj-dcgQk I I 

1 < Vk/" ^ K < l1 



< 



j=i v T ' 0<i<i fc 

k-j+degl 
Cl ( n _ -|- ^dffe-i+dogQj-dcgQfc) 

^ j (n - jfc + 

jeB v ^ 7 



(n- k+ l) d 



5>n 



fe-j+dog Qj-degQ fe , 

\ - \ak,: 



jec y ' a<i<i k " 

where the last three sums tend to zero when n — > oo by the above arguments 
(and the last one due to Theorem 1). 

Thus for all sufficiently large n there exists a positive constant c' such that 

where c' — ► 1 when n — > oo. If the set A contains precisely one element, then the 
sum in (|15p consists of one single term, and we are done: there exists a positive 
constant c = ( c y2^ 1/(fe ~ i+desQ;/ ~ desQ,t) such that r n > c a (n ~ k + l) d for all 
sufficiently large n, and thus limn—^ inf (r n /n ) > c$. 

But clearly, for some operators A will contain more than one element. If this 
is the case define m := mirij g A(fc — j + degQj — degQfe) and denote by j m the 
corresponding j for which this minimum is attained. Using the upper bound 



31 



T n < ci(n — fc + l) d we then get the following ine qualitjH from ([15]): 

/ „ \ i'-j+dcgQj-degQfc 



v (n - £: + l) d 



if, 



3m \ (n-k + l) d 



E *j 



(n-k + l) d J \(n~k+l) d 



fc-j+dcg Q 3 -degQj.-m 



< 



i£A\{j m } 

r » \ \ ~> y / r " \ fc-j+dcg Q, -deg Q 



(n - k + l) d 



(n-k+ l) d 



K 3m + E K r c " 

jeA\{j m } 



k-j+deg Qj -deg Q k -m 



■K 



where K > 0. Thus r„ > (t^) 1 ^™^ — fc + l) d for all sufficiently large n, and 
therefore there exists a positive constant cq — (l/K) 1 ^" 1 (recall that c' 1 
when n — > oo) such that linin^oo inf £f > cq. □ 



9 Consider for example the operator T = zD + D 2 + zD 3 + 2;_D 4 . Here by Lemma 3: 
2 4-j 3-j+dcgQ 3 3 r r 



(n-3) 4 -J (n-3) 3 (n - 3) 2 (n - 3) ' 

where r n is the largest modulus of all roots of the unique and monic eigenpolynomial of T. For 
this operator d = 1 and we see that (4 — j)/(3 — j + deg Qj ) = d for j ' = 1 and for j = 3. Now, 
assuming that r„ < ci(n — 3) holds for some positive constant ci for large n, our inequality 
becomes 



1 < 



(n - 3) 3 (n - 3) 2 (n - 3) 



r„ cf(n-3) 2 ci(n-3) r. 



(n-3) (n-3) 2 (n-3) 2 (n-3) 

= ( 8c 2 + 2 )^ + ^ 
(n — 3) (n — 3) 

where the last term tends to zero as n — > oo. Thus r n > co(n — 3) for sufficiently large choices 
on n, where cq = l/(8c 2 + 2), and hence limn—^oo inf(r n /n) > cq. 
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6 Open Problems 



1. The main challenge is to obtain a complete proof of the Main Conjecture, 
see Introduction. This proof requires both the sharp upper and lower bounds 
of the largest root. The upper bound can apparently be obtained by a detailed 
study of the corresponding Riccati equation at oo. If the Main Conjecture is 
settled then to achieve its corollary (the Cauchy transform equation) one can 
use a technique similar to that of [5] to prove the basic assumption, see Section 3. 

2. As suggested by one of the referees, estimates similar to that of Main 
Conjecture can be formulated for the sequence of roots 2 n ,i of p n such that 
lim^oo = a where < a < 1 and \z n<i \ = r n ,j. 

3. In the case of orthogonal polynomials there is a number of results describing 
the growth of the largest modulus r n of the roots as an expansion in powers of n, 
see e.g. [llj and [12] and references therein. In the present paper we conjectured 
the form of the leading term of r n in our more general setup. As suggested by 
one of the referees, the question about the lower terms in the expansion of r n 
is natural in our context as well. 

4. Operators of the type we consider occur in the theory of Bochncr-Krall 
orthogonal systems, i.e. families of polynomials which are both eigenfunctions 
of some finite order differential operator and orthogonal with respect to some 
suitable inner product. A lot is known about the asymptotic zero distribution 
of orthogonal polynomials, and by comparing such known results with results 
on the asymptotic zero distribution of eigenpolynomials of degenerate exactly- 
solvable operators, we believe it will be possible to gain new insight into the 
nature of BKS. We have previously used our results from [2] to prove a special 
case of a general conjecture describing the leading terms of all Bochner-Krall 
operators, see [3]. Another problem relevant for BKS is to describe all exactly- 
solvable operators whose eigenpolynomials have real roots only. 

5. Numerical evidence indicates that the roots of the scaled eigenpolynomi- 
als fill certain curves in the complex plane. The support of the limiting root 
measure fj, seems to be a tree. This is the case for the non-degenerate exactly- 
solvable operators which we treated in [2] , but then without such a scaling of the 
eigenpolynomials. By a tree we mean a connected compact subset T of C which 
consists of a finite union of analytis curves and where C \ T is simply connected. 
The (conjectural) algebraic equation satisfied by the Cauchy transform contains 
a lot of information about /i, and it remains to describe its support explicitly. 

6. Conjecturally the support of the asymptotic zero distribution of the scaled 
eigenpolynomial q n is the union of a finite number of analytic curves in the com- 
plex plane which we denote by St, i.e. 2<r = supp ji, where \i is the limiting 
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root measure of q n . Then the following conjecture seems to be quite plausibleB 

Conjecture 1. [Interlacing property] For any family {q n } of appropriately 
scaled eigenpolynomials of a degenerate exactly- solvable operator, the zeros of 
any two consecutive polynomials q n +i and q n interlace along St for all suffi- 
ciently large integers n. 

When defining the interlacing property some caution is required since the 
zeros of q n do not lie exactly on St . Thus identify some sufficiently small neigh- 
bourhood N(St) of St with the normal bundle to St by equipping ./V(St) with 
the projection onto St along the fibres which are small curvilinear segments or- 
thogonal to St- We then say that two sets of points in ./V(St) interlace if their 
orthogonal projections on St interlace in the usual sense. If St has singularities 
one should first remove some sufficiently small neighbourhoods of these singu- 
larities and then proceed as above on the remaining part of St- Conjecture 
1 thus states that for any sufficiently small neighbourhood N(Et) of St there 
exists a number uq such that the interlacing property holds for the zeros of q n 
and q n +i for all n > n . We conclude this section by showing some pictures 
illustrating the interlacing property. Below, small dots represent the roots of 
q n +i and large dots represent the roots of q n for some fixed n. 



-0.2 
-0 . 4 



-0.6 -0.4 -0.2 0.2 0.4 0.6 

T = z 2 D 2 + z 3 D 3 + zD 5 
roots of <?25 and 524- 



-1 -0.75-0.5-0.25 0.25 0.5 

T = zD + z 2 D 2 + D 3 
roots of (/20 and qi$. 



-0.25 
-0.5 
-0. 75 



-1.75-1.5-1.25 -1 -0.75-0.5-0.25 



-0.8 -0.6 -0.4-0.2 0.2 0.4 



T = zD + zD 2 + zB z + zD 4 + zD 5 , T = z 3 D 3 + z 2 D 5 + zD 6 , 
roots of 523 and 522- roots of (720 and qig. 



°The question concerning interlacing was raised by B. Shapiro. Also see pQ. 
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